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Cluster Dynamical Mean Field Theories are analyzed in terms of their semiclassical limit and 
their causality properties, and a translation invariant formulation of the cellular dynamical mean 
field theory, PCDMFT, is presented. The semiclassical limit of the cluster methods is analyzed by 
applying them to the Falikov-Kimball model in the limit of infinite Hubbard interaction U where they 
map to different classical cluster schemes for the Ising model. Furthermore the Cutkosky-t'Hooft- 
Veltman cutting equations are generalized and derived for non translation invariant systems using 
the Schwinger-Keldysh formalism. This provides a general setting to discuss causality properties 
of cluster methods. To illustrate the method, we prove that PCDMFT is causal while the nested 
cluster schemes (NCS) in general and the pair scheme in particular are not. Constraints on further 
extension of these schemes are discussed. 



I. INTRODUCTION 

Dynamical Mean Field Theory (for a review see 0) has been very successful in providing a non perturbative 
approach to strongly correlated Fermi systems. It describes both the localized and the itinerant limit and has yielded 
non perturbative insights into the finite temperature Mott transition 1]. This approach has been combined with 
realistic electronic structure methods such as LDA and GW and has been applied successfully to numerous materials 
PI- In spite of these successes, several limitations of the single site I^MFT approach are now apparent. For example, the 
self-energy is k independent by construction, so the method cannot describe independent variations of the quasiparticle 
residue, the quasiparticle lifetime and the effective mass. Furthermore, the single site nature of the method precludes 
it from treating more exotic orders with order parameters involving several sites, such as dimerization, staggered flux 
or d-density wave, and d-wave superconductivity. 

To overcome these limitations, various extensions of DMFT have been proposed. For a disordered system, one can 
set up a functional integral formulation, with DMFT as a saddle point, leading to a natural loop expansion l3|. A 
formulation of these ideas for a clean system is still lacking. A different extension, is the Extended DMFT in 
which the DMFT Ansatz is applied simultaneously to bosonic and fcrmionic degrees of freedom. This approach does 
not describe a fc-dependent self-energy but incorporate diagrams involving longer range interactions into the DMFT 
equations. 

A different idea is based on truncations of the Baym Kadanoff functional. In the full theory, this is a functional of 
the full Green functions. DMFT is obtained by restricting it to local Green functions only, setting the non-local Green 
functions equal to zero. A natural extension is to restrict the functional to local and nearest neighbor Green functions. 
This pair scheme was introduced independently by Ingersent and Schiller j5|, and by Georges and Kotliarjlj. The 
latter authors observed non causal behavior in an Iterated Perturbation TheoryJIPT) solution of the Hubbard model, 
but the origin of this problem was not elucidated. Zarand et al. suggested 6] that the pair cluster approach was 
causal, and that the difficulties encountered in the solution were related to the impurity solver rather than to the 
scheme itself, but no conclusive proof of this statement was presented. This pair method requires a simultaneous 
solution of multiple site impurity problems. 

A different direction was pursued by Jarrell and collaborators with the introduction of Dynamical Cluster Approx- 
imation [DC A ) 0, whose main idea is to discretize uniformly momentum space. This approach was shown to be 
manifestly causal. It can also be formulated in real space ^] (See also Appendix^. A different approach, the cellular 
DMFT or CDMFT, motivated by applications to electronic structure, was introduced in In this approach, the 
many body problem is truncated by introducing a finite basis set of orbitals to truncate the self energy. It introduces 
the cluster self energy and the lattice self energy as independent entities. This method was tested in a soluble model 
0, in the Falikov-Kimball model and in the one dimensional Hubbard model 0. These papers developed 
the CDMFT ideas from a real space perspective. In this paper we develop this approach from a momentum space 
perspective, to maintain periodicity in the self consistency equation. The importance of including this periodicity 
was underlined by Lichtenstein arguing that the lack of periodicity of CDMFT could unfavor phases, like 

the d-wave superconductivity, in which the order parameter lives on links. Even though this point has not yet been 
elucidated, it is certainly desirable to have a generalization of CMDFT that is translation invariant and causal in its 
use of the lattice self-energy in the self consistency condition. 
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Cluster DMFT methods have not yet reached the level of understanding of their single site counterparts. While the 
single site DMFT can be unambiguously formulated, cluster DMFT methods are more diverse and therefore require 
more detailed methodological investigation, since the virtues and limitations of the various cluster schemes are not 
apparent yet. This paper is a contribution in this direction. 

There are several important principles that a cluster method should satisfy: i) Given that cluster _DMFT approxi- 
mations are intrinsically basis dependent, a cluster method should be formulated in a general basis set. This flexibility 
is important, because for a given problem, one could carry out the cluster DMFT study in the basis which is most 
suitable for the system in question, ii) It should have an effective action formulation, namely it should target the 
calculation of a specific correlator function. Hi) It should yield causal Green functions, iv) It should be able to 
capture the various order, including those which break the translation invariance. v) It should converge rapidly as a 
function of the discretization parameter for the observable that one is targeting. It is possible, that different cluster 
schemes may converge better for different observables. A better understanding of these elements of a cluster method 
is desirable. 

In sectionini we present our translation invariant generalization of CDMFT where the lattice self energy participates 
in the self consistent equation. This is a new cluster scheme, and we formulate it in a way which allows a comparison 
with CD MF T a.nd DCA . We also discuss in section III CI a different class of schemes, the nested cluster schemes^ which 
require the simultaneous consistent solutions of impurity models of different size. 

In section Unl in an attempt to clarify the nature of the various cluster approximations, we study their classical limit 
in the Ising limit of the Falikov-Kimball model. The various schemes then reduce to classical cluster approximations 
to the Ising model. This analysis elucidates their physical content at the classical level and it allows a simple and 
clear comparison. This study complements other comparisons of cluster schemes against more exact treatments for 
specific models 0, ^ . The classical limit of EDMFT was discussed in Q . 

In section IIVI we analyze the issue of causality of cluster methods from a diagrammatic perspective using the 
Cutkovsky-t'Hooft-Veltman rules. We rederive them for non translation invariant cases using the Schwinger-Keldysh 
formalism. We end up with a general setting to analyze the causality of cluster schemes. Within this framework we 
show that the pair scheme is not causal, and elucidate the origin of the problem. Furthermore we provide a simple 
proof of the causality of CDMFT and DCA (which had already been proved by other methods) and justify the 
causality of the periodi c ge neralization of CDMFT. Finally we make a relationship with the early paper of Mills et al. 
on disordered systems ilSj: in fact the coherent potential approximation (CPA) is a particular case of a DMFT. The 
origin of violations of causality, and their possible cure for the generalizations of CPA were clarified by Mills l^. 

II. DESCRIPTION OF THE CLUSTER METHODS 

A fairly general model of strongly correlated electrons contains hopping and interaction terms. It is defined on a 
lattice of M sites in d dimensions, and we divide the lattice in (M /L)'^ cubic clusters of Sc = L'^ sites (more general 
forms can also be considered). We denote with the internal cluster position and with R„ the cluster position in 
the lattice (therefore the position of the i-th site of the n-th cluster is i?„ + r^). The lattice Hamiltonian is expressed 
in terms of fermionic operators ^ and CR^_p and can be written as: 

^ ia,/3(i?n - i?m)4„,aCi?,„,/3 + ^1 t^",/3,a',/3' ({^})4„ ,/34^^, C/J^, ,/3' (1) 

n,a,rn,l3 n,a,m,l3,n' ,a' ,m' 

where a = i,a and a is an internal degree of freedom (i.e. a spin, spin-orbital or band index). Note that in all this 
paper, we will denote the position and momentum on the (original) lattice with lowercase letters (r and k) and the 
position and momentum of the superlattice with uppercase letters (i? and K). Different dynamical cluster methods 
for strongly correlated have been introduced in order to obtain an approximate solution of able to capture the 
effect of short range (dynamical) correlation and to describe the self-energy k-dependence. 

In the following, we introduce a new scheme, the Periodized Cluster Dynamical Field Theory {PCDMFT), formu- 
late CDMFT and DCA in real and k space and the nested cluster schemes (NCS), i.e. the pair scheme and its 
generalisations. CDMFT is a real-space cluster : the lattice is divided in a superlattice of cells and the scheme is 
basically the ZJMFT equations on the superlattice P|. On the other hand, DCA is a reciprocal space cluster, where 
the self-energy in momentum space is approximated by step function around a few points, which we will denote by Kc 
and are identified to the momentum of the (periodic) cluster. We will show that PCDMFT is a natural generalization 
of both schemes, from a real space and a fc-space perspective. 
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A. Real space perspective 



Cluster methods for strongly correlated electrons {CDMFT, DC A , PCDMFT) can be divided in two steps. (Nested 
cluster schemes will be considered separately). 

The first step is the computation of the local cluster propagator Gca^viT) and the cluster self-energy Tp from an 
effective action containing a Weiss dynamical field Gg_^^j^(r, r') and the intra-cluster interaction : 

S,s^- drdr'cl^ir)Go,l,,Ar, t')c„,{t') + / drUo^p-^siR = 0){cicpc\cs){T) (2a) 
J Jo Jo 

Gc.^At) = - {Tc,,iT)cUO))s.„ (2b) 



- Gq-i - G-i (2c) 



where fj,, v are indices for cluster sites. In the following, we will often concentrate on square clusters of linear size L 
on a d-dimensional square lattice, although many of the results can be generalised easily. Hence, /i, v will also denote 
the position of the cluster sites on the lattice : fj.ji' G C — {0, . . . ,L — 1}'^ where d is the dimension (intersite distance 
is normalised to one). We will denote by Sc the cluster's size and by C the set of cluster points; a is the spin index 
as mentioned earlier and a, (3, 7, 5 are double indices gathering the spin and the cluster index. 

The second step consists in recomputing the Weiss field using the value of the self-energy obtained by the first step 
and then iterating until convergence is reached. The real difference between cluster schemes is how the second step is 
performed. 

CDMFT is a direct generalization of DMFT to cluster in real space and consists simply in rewriting the DMFT 
equations in a matrix form (/i, v indices are omitted) taking as elementary degrees of freedom all the cluster fermionic 
degrees of freedom: 

\KeitB.z. + ^ ~ H-^) ~ ^ («^") / 

where t^^{K) is the hopping expressed in the superlattice notations, with K in the Reduced Brillouin Zone (R.B.Z.) 
of the superlattice. (See Eauation lAll in Appendix). Note that from now on the sum over K means always the 
normalized sum. When convergence is reached and the cluster self-energy has been obtained, the translation invariant 
lattice self-energy E'***' {{hj) denotes a site on the original lattice) is computed by the formula: 



^latt 



(c.) = i- J2 ^U^) (4) 



Other estimators (that still preserve causality) have been introduced in Note that in CDMFT, the cluster and 
the lattice self-energy are two different quantities (2j : one has to first solve the cluster problem, and only at the end 
to compute the lattice self-energy E^*^**. In CDMFT, E^^" does not enter the self-consistency condition. 

DC A is more naturally formulated in Fourier space 0. However it can also be recovered from a real space 
perspective @ just changing the hopping matrix in eq. Q to (See also Appendix IX)) : 

i^!;^{K) = t^,{K) expi-iK ■ (^ - i^)) 

In this case if the self-energy is a cyclic matrix (translation invariant within the cluster) then eq. reduces to 
(upon diagonalisation) : 

G^\K,,iu;^) ^ ( V . ■ ^ . , ^ +EC(X„za;„) (5) 



\KeR.B.Z. 



lOJ, 



that is the standard DCA equation introduced in 0. Note that if the self-energy is cyclic then the Weiss field 
computed by Q will be cyclic, so this property is preserved within the self-consistent loop. In DCA , contrary to 
CDMFT, there is no distinction between the cluster and the lattice self-energy. Because it is formulated in fc-space, 
DCA is also naturally translation invariant. 

We are now ready to define the new scheme, PCDMFT for Periodized CDMFT. The simplest definition is to take 
j^iatt jj-^ place of E^ in the self-consistency condition. This is very close to the scheme proposed by Lichtenstein and 
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Katsnelson in with the big difference that PCDMFT is causal, as will be proved below. Thus the equation 
relating the self-energy to the Weiss field are: 



T,iatt{k, iujn) = 4- X! exp(-iA: • - i^)) (6a) 





1 






= E 




k 




= G- 



^ — ik-fi^ik-v 

(6b) 



iuJn + t{k) - i;ia"(A;, iuJn) 

' + (6c) 

where k is in the Brillouin zone of the original lattice. 

The three schemes (CDMFT, DC A , PCDMFT) can be summarized into the same matrix equations : Eq. ||2Jl and 



-1" 



-1 



Go\iLUr,)^i J2 iiiUn + f^-isiK)~j:s{K,iu;n)] +S^(ic^n) (7) 

XKeR.B.Z. ^ ^ / 

where the difference between the three schemes is enclosed in the value of ^5 and of E5 that enter in the self-consistency 
condition. 

Let us now turn to translation invariance breaking phases, where translation invariance is conserved on the su- 
perlattice (e.g. antiferromagnet, charge density, "stripes" if the cluster is big enough). CDMFT can describe such 
an order by construction since it does not require the translation invariance. When solved numerically (e.g. with 
Quantum Monte Carlo method), the translation invariant solutions are often found to be instable towards the ordered 
one (we have explicitely encountered this phenomenon for example for AF and charge density wave). 

DCA and PCDMFT require translation invariance, therefore they need to be generalized to handle such an order. 
For DCA there are two solutions for this problem, for example for antiferromagnetic order : i) keep a reciprocal s pac e 
formulation with a Reduced Brillouin Zone and introduce some correlation between k and k + Q, with Q = (tt, tt) 16j; 
ii) use the real space formulation introduced in ,^J, where translation invariance in the cluster can be broken, and look 
for an antiferromagnetic solution. It is shown in Appendix ^ that the two approaches are equivalent. However, i) 
requires to anticipate the appearance of ordered phase, i.e. to adapt the cluster scheme for the order to be described 
: one needs to use a special setup for antiferromagnetic order, another for a more complicated order, whereas ii) does 
not require to anticipate the order : it will show up automatically solving the real space DCA equations with no 
need to generalize the scheme (provided that the cluster is big enough to contain the unit cell). The same numerical 
code will produce a translation invariant solution, or an antiferromagnetic one, a stripe-like one. In particular, since 
translation invariant solution are often found to often be numerically instable, one can not miss a ordered phase with 
this approach. Therefore, the real space formulation of DCA is the best solution from a practical point of view. 

The PCDMFT is unfortunately more complicated to generalize. In the following, we will focus on square clusters 
on the square lattice. It is useful to introduce a slightly more general formula for E''^** for a bipartite lattice. The 
lattice self-energy is a sum of the cluster self-energy put at all possible positions on the lattice. We can rewrite Q in 
a more transparent way, as a sum over all possible shifts of the cluster : 

5e{0,...,L-l}<i 

where \_x\ is the integer defined by \_x\ < x < [xj -I- 1 (for a vector it has to be understood component by component), 
d is the dimension, (5 is a d-dimensional shift vector, the bar denotes the modulo reduction by L component by 
component (adding (5 is a circular shift in the cluster). 

The idea of the generalization is then simple. In a phase that breaks translation invariance (we take antiferromagnet 
as an example), we have multiple solutions (two for AF), denoted by an index a (a = 1, 2 for AF). In order to respect 
the order, we need to compensate the shift in the cluster by the change of solution : in the AF example, use one 
solution on one sublattice, another on the other sublattice. The formula is then 

^ Se{o....,L-i}'' 

^ ' \2 if Eto'5.-l [2] ^ ' 

For the AF phase and when the original model is SU(2) invariant, the solution a = 2 is simply obtain from a = 1 
by a spin-flip. This does not need to be true in general, either for more complicated order or the Falikov-Kimball 
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model that we will use in Section UTTl Let us note however that this generalization is more complex than for the DC A 
scheme from a practical perspective : the expression of S'^** strongly depends on the order to be described and on 
the form of the cluster. 

B. Reciprocal space perspective. $ derivation. 

We now explore the relation between PCDMFT, CDMFT and DC A from the reciprocal space perspective. In order 
to do this, we use the generating functional formulation of DMFT IJ. Let us recall that one can define a functional 
r(G): 

r(G) =TrlogG-TrG^iG + $(G) (10) 

where Gq = [i'jJn — t + fi)^^ is the bare propagator and <i>(G), the Baym-Kadanoff functional, is the sum of all 
the vacuum two-particle irreducible diagrams constructed with the propagator G and the interaction vertices. The 
solution of the stationarity equation of (|10|l is the real propagator of the full interacting theory. DMFT, as well as its 
cluster generalizations, can be seen as an approximation onto $(G). Indeed one obtains the iJMFr approximation 
restricting ^(G) only to single site propagator (equating to zero all the non single site propagator), i.e. ^dmft{G) = 
'^^^{Gii\Gij — 0) or, equivalently, neglecting the momentum conservation at the vertices. These two procedures, 
that are equivalent at the single site level, represent two different routes in order to obtain cluster generalization of 
DMFT. 

CDMFT can be obtained by a natural real space extension of the DMFT approximation on 

^CDMFt{G) = HG^M;uM\GpM;XRf = 0) (11) 
R 

It's easy to show that the stationarity equation corresponding to this choice of $ gives back the CDMFT equations. 
Note that the summation of the infinite series of diagram of ^cdmft{G) is performed by the cluster impurity solver 
like in the DMFT case. 

DCA is formulated more naturally as an approximation in Fourier space: instead of neglecting completely the 
momentum vertex conservation one puts a coarse-grained delta inside the complete diagrammatic series. 

<S>dca{G) = N,s '^{Gimu^k^^k^^k^^.^^^UocAik.MMM (12) 

where C/DCA(fci, fe, fcs, ^4) = 5K^(ki)+Ka(k2),Kc{k3}+Kaiki)/Nss, Nss is the number of clusters (number of sites divided 
number of sites per cluster) and Kc{k) is a function that for each k gives the center of the j^]'^ cube to which k 
belongs 0|- 

PCDMFT has been introduced previously from the real space point of view like a natural generalization of CDMFT 
in which one puts the lattice self-energy inside the self-consistent loop (still preserving the causality properties of 
CDMFT). In the following we shall show that, from the functional point of view, PCDMFT can also be seen at the 
same time as a translation invariant formulation of CDMFT and as a generalization of DCA . Up to now we restricted 
the discussion to the standard completely localized basis set for simplicity. However, CDMFT and PCDMFT can be 
formulated in a general basis set. This flexibility is important, because for a given problem, one could carry out the 
analysis in the basis which is most suitable for the system in question. As a consequence in the following we will 
derive PCDMFT for a general basis. 

We shall show that there are two different procedures leading to the same formulation of PCDMFT. Let us focus 
first on the one which show that PCDMFT is the natural translation invariant formulation of CDMFT. Call (paa the 
basis function used to define the cellular DMFT and (pua their Fourier transform. Note that they are normalized in 
such a way that J2i \'PRa{i)\'^ = Ij X]fc \'?Ra{k)\'^ = 1. The particular case of the standard completely locahzed basis set 
corresponds to ^Ra{k) = exp{—ik-xiia)/VN {N is the total number of sites). CDMFT in a general basis is obtained: 
i) keeping only the intra-cluster interaction: U ii^a,R2[3 .,Rzi ,RiS — > URa,R(3,R-j,RS and ii) making the approximation Hll|) 
discussed above to the Baym-Kadanoff functional. 

PCDMFT cam be obtained making the approximations i), ii) and imposing the translational invariance of the original 
problem. This is performed expressing the propagator in the new basis set in terms of the translational invariant 
propagator in the original basis set: 



GRa,R'f3{i0) = Y(pB,a{k)*G{k,Uj)(pR'p{k) 
k 



(13) 
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Plugging this expression inside the functional leading to CDMFT we get the PCDMFT functional: 

T{G{k)) = Trlog(G) + Tr((iL^„ - t{k) + iJ,)G) + Nss<^pcdmft{U , G) (14) 

where $ is the Baym KadanofF functional expressed in terms of G{k) and obtained replacing the original matrix of 
interactions U(ki, A;2, ^3, ^4) with 

U{ki,k2,k3,k4) = ^ t/0a,0/3,07,05'?0a(fcl)^S/3(fc2)^07(^3)^S5(^4) (15) 

where we made use explicitly of the translation invariance to put i? = 0. 

Extremizing this functional with respect to G{k) gives the PCDMFT equations together with the form of the lattice 
self energy: 

GRa.^R/niuj) = V ^*Rc,ik) 1; 77 — ^VRfs{k) (16a) 

^ LO — tk — 2^latt\K,Uj) 

^iattik,uj) = ^*j^^{k)lllf^i:pRp{k) (16b) 

afJR 

where is the self-energy obtained from a cluster impurity problem characterized by a propagator ^ = Gr^^rp 
and the interaction matrix Uap-yS- 

As discussed before there is another procedure to get the functional formulation of PCDMFT encoded in (|14|l . As 
the DMFT approximation can be obtained just by neglecting the momentum conservation at the vertex, one can 
interpret (|15|l as improvement to the DMFT approximation in which the vertex is replaced by U{ki, k2, k^, k^) = 

{kl)'Pop{k2)'POi{k3)ipos{k4) where Ua,f3,j,5 = C^0a,0/3,07,05- 

At this point a natural question is why we have chosen this particular U and, more interestingly, could other choices 
lead to a better approximation and what is the procedure to find the "best" {/? Unfortunately we don't have clear 
answers to these questions. A partial answer to the first one is that the choice of U leading to PCDMFT is the one 
that we obtain applying the least square minimization to U{ki, k2, k^, k^) — U{ki, k2, fca, k^). So in a certain sense it 
provides the best approximation to U within the chosen basis set. It's however important to notice that is far from 
clear that this is a good criterion to select the "best" tj . With respect to these remarks is particularly interesting 
to note that DCA can be obtained taking Ua,f3,'y,5 = t^oa,OQ,Oa,OQ^a,/3<57, (5^(5,7 (in the case of a complete diagonal 
interaction on the original lattice) and (pRi3{k) = exp{—iKc{k) ■ XRj3)/^/N 0. 

This suggests interesting interpolation between DCA and PCDMFT. A particularly simple one consists in dividing 
the Brillouin zone in l'^ squares (like in DCA for a cluster of linear size I) and taking ipR^{k) — exp{—iFi{k) ■ x Rf^) / \/N 
where Fi{k) is a function such that for each k gives the center of the [j, j]'^ cube to which k belongs (note that 
fi = 1, . . . , Lc). When / — Lf. one gets DCA , whereas when / = oo one gets PCDMFT. 

Let us finally note that this functional derivation of PCDMFT is crucial to prove simply its causality properties. 
Indeed this scheme is causal because it is a composition of two causal steps. Clearly starting from a causal cluster 
self-energy one obtains a causal cluster propagator through (|16all6b|l . Furthermore we will show in Section Hvl that 
plugging the causal cluster propagator inside the diagrammatic series corresponding to the PCDMFT functional gives 
a causal cluster self-energy. 



C. Nested cluster schemes 



Another natural generalization of DMFT to the cluster setting is to apply the ideas connected to the cluster 
variation method (CVM) of classical statistical mechanics to the Baym Kadanoff functional T?, Ts'l. These schemes 
are also natural generalizations of the n site CPA to interacting electrons. The approach is defined by selecting 
a set of maximal (namely they are not included in each other) clusters of sites. We denote by F the set of maximal 
clusters together with all its subclusters, and by the restriction of the Baym KadanofF functional to Ga{i,j) , with 
Ga{hj) = G{i,j) if i and j belong to a and Ga{i,j) = otherwise. $q, is defined recursively in terms of the 
by $a = X]/3Cq which can be inverted by the Moebius formula = S/3Ca(~-'^)"'^~"°^/3 "^'^^^ '^a the number of 
sites of cluster a. An approximation scheme is uniquely fixed once a set of maximal clusters is chosen. If the chosen 
set is invariant under translations, e.g. the set of all plaquettes, we construct a cluster scheme which is manifestly 
translation invariant by truncating the full Baym Kadanoff functional 



OLfT 



(17) 
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Differentiation of $ yields a translational invariant self energy, which requires the solution of several impurity 
problems. The subclusters of a maximal cluster are generally related among each other by the operations of the 
crystal group, and fall into different equivalence classes. To compute the lattice self energy one needs to solve several 
impurity models. One for each representative of inequivalent subclusters of the maximal cluster. 

When we take as a set of maximal subclusters the set of all nearest neighbor pairs of the lattice we obtain the pair 
scheme (or two impurity scheme) 0, 0] which we discuss in more detail for completeness. It is convenient to go back 
to a lattice notation where lattice sites are denoted by i and j and not by The approximation to the $ functional 
can be written in terms of $i and $2 which are the Baym Kadanoff functionals of a one and a two impurity problem. 

<^pa^r = (1 - z) ^ $i[G,,] + ^ $2 [G,« , Gjj , ] (18) 
i <ij> 

Differentiating this functional gives an equation for the local self energy and the nearest neighbor self energy. 

+ z — — — (19a) 



6Gu SGii \ SGii 6G 



5Gij 



(19b) 



The diagrammatic interpretation of the first equation is transparent : each diagram that involves a link is properly 
counted in the solution of a two impurity model. The diagrams for the local self energy at site i which involve only 
the Green function at site i, can be obtained from a one impurity model with the correct counting. The contribution 
of the diagrams for the local self energy which involve the Green functions of a pair of sites, can be obtained by 
subtracting the contributions from the one and the two impurity self energy and multiplying by z, which results in 
equation p9a|l . From Sjoc and we can construct the lattice self energy. 

E'"^" = ^lociiuJn) + t{k)^nn{i^n) (20) 

and close the equations by requiring the self consistency condition imposed by the Dyson equation 0. 

i^n - t{k) - J^lociiuJn) - t{k)Y.nn{iuJn) ^ 

( G G \ 

This can be expressed in the matrix notation, with 2x2 matrices G — \ I and K in the reduced Brillouin 

zone : 

G= J2 {i^n-t{K)-Y}-'\K)y^ (22) 

KeR.B.Z. 

The generalisation of this scheme to antiferromagnetic order is presented in Appendix^ An important feature of this 
scheme is that the Green function of the one site problem coincide with the diagonal part of the Green function of the 
two site problem. This is a general "nested" structure for these schemes, hence the name "nested cluster schemes" . 
We will see in section Hill that this property leads to a quantitatively good classical cluster scheme. 

D. Hartree-Fock terms 

Another important issue is the treatment of longer range interactions within CDMFT, and PCDMFT. In this 
context it is worth noticing that the Hartree Fock contribution to the Baym Kadanoff functional (|10|l 

'^hf[G] ~ ^ Ua[i^s{Ri, R^, R4){G pa{R2Ri)Gs^{RiR'i) — Gsa{RiRi)Gfj^{R2R3)) (23) 
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induces a self energy which is frequency independent and therefore does not cause problems with causality and can 
be evaluated with little computational cost. So it is convenient to separate $ = ^hf + ^dyn, and apply the cluster 
DMFT truncation only to ^dyn, and to the self energy it generates while treating the Hartree contributions exactly. 
More precisely, one can treat with Hartree-Fock terms that connect the cluster to the exterior only, to avoid a double 
counting problem. 

This observation is particularly relevant in the treatment of broken symmetries induced by non local interactions as 
exemplified in the study of the transition to a charge density wave in the extended Hubbard model in one dimension 
which was studied in ref Hill. 



III. CLASSICAL LIMIT 



In this section, in an attempt to clarify the nature of the various cluster approximations, we investigate analytically 
the large-C/ limit of the Falikov-Kimball model, which reduces to the classical Ising model in that case. The Falikov- 
Kimbal model is defined by the Hamiltonian : 

H^J2 *y-clc»<r + U (n,| - i) - i) (24) 

where denotes nearest neighbors, = t and Uji^ ~ 0. We consider the particle-hole symmetric case (/i 0). 
This model has been studied a lot for its own interest (for a review, see 20]), but we will use here as a tool to derive 
a classical limit for the various cluster methods. This completes the QMC study of this model with CMDFT and 
DCA 10]. Indeed, in the limit [/ — > oo with /? oo, (3 = pt/U fixed, it reduces on the lattice to the Ising model at 
temperature : 

H=Y^ jf^S^S] (25a) 

Jl^f = ty/2 (25b) 

with Sf = (cJ|Ci| — c^cu) (Ising spins). Note that a factor t/U has been absorbed in the definition of p. The proof 
is analogous to the standard reduction of the Hubbard model to the Heisenberg model : since down electrons are 
quenched, there is no possible exchange between two spins at different sites, which implies that the interaction is 
Ising-like. We will now take the classical limit of the various cluster schemes, using their common expression ^ 
and O in order to obtain classical cluster approximations of the Ising model. For the reader not interested in the 
derivation, we summarize our findings in section Fill Bl 



A. Derivation of the classical limit 



We first present the derivation of the classical limit of CDMFT, DCA and PCDMFT. Our results on the classical 
limit of the Pair scheme are presented but derived in appendix lUl 

In the following, we will concentrate ourselves on square clusters on the square lattice. For clusters with an odd size 
it is necessary to slightly generalize the CDMFT equations to get antiferromagnetism in the usual way, distinguishing 
two sublattices To avoid cumbersome notations, we will discuss that point later. The key point is of course that 
tiji = makes the model partially solvable in the various schemes. 

A priori, we have an effective action given by 0. Since t| = 0, we can show self-consistently that G04 is diagonal. 
Indeed, if Gqi is diagonal, because of the form of the n-^ni vertex, Gi is diagonal (at all order in U), and so is Sj. 
In all schemes, Sj^** is then diagonal and independant of K. Therefore from we get for CDMFT, DCA and 
PCDMFT. 

(Gox) \iw„) = iujn + S| (iwn) - Sj^"(ia;„) = iu-a 

(For the pair scheme, see Appendix IdI) . Therefore we can consider rt^ as a classical variable and compute the Green 
function for the up electrons, solving the effective action : 

^eff - - / f dTdT'cl^{T)G^l^,{T,T')c,^{T') + f dTcl^{T)drC^^{T) + [/(tI^tW " \) ("M W " \) (26) 
J J Q J Q 
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where /j,, v are cluster indices. For fixed n^, the action for the up electrons is Gaussian, which leads to : 



G,.= E ^^^Y^iG-.l-{n,,-\)U5,A (27a) 
Z{{n,,}) EE exp(^Trln(^Go-:, - (n^^ - i)C/<5^.) j e^^/^^.""^ (27b) 

{npi=0,l} 

The computation of the large-[/ limit is organized in two steps : first, we find the expansion of Gq, and second 
we show that in this limit the effective action Ij27bp becomes the action for a classical Ising cluster with mean-field- 
like terms. Moreover, the study of the large- ?7 limit requires an expansion in the limit U ^ oo and oj ^ oo, with 
X = 2uJn/U fixed. Indeed there is a strong lo dependence at the scale as seen already in the atomic limit. We will 
see that the final result is not determined by the a; — > limit. 

First, we use an Ansatz for Gq, that we will prove to be consistent 

GoHiuJn) ^ i^x - Aix) (28) 

with A of the order 1 in 1/U (plus subdominant terms). Using the expansion 

(A - A)"^ = A-^ + A-^AA-^ + ■■■ 

with A^jy — {Sfj, + ix) Sfj_ = 1 — 2n^|, we expand G| in H27a|l and obtain from the following expansion for 

the cluster self-energy of the up electrons S*-" : 

S^.-S^;r + <5I]^. (29a) 
^^T-^{j^-)^Oil) (29b) 

c =(l-5 )A ^"^^ ■ - ^ ^ 

-r ix 

where j]^'^''*^ (resp. <5S^^) is the diagonal (resp. off-diagonal) part of the matrix, the brackets denote the average 
over n^i or with the weights defined in (|27a|l . The averages and correlations of are computed using S*^ = ±1 
and solving for the probability of the spin to be ±1 as a function of the correlations: 



-1 



^I^;r.-(l-MA..^)^+0^-) (29c 



{K) = ^-^A_^ (30a) 



, > {S,S,)^^x{{S,) + {S,))-x^ 

{h^K) = f -2 (30b) 

(1 + x^j 

For (/i, v) nearest neighbors in an antiferromagnetic phase, as well as in the paramagnetic phase, we have (5*^) + {Sv) — 
and thus: 

jh^K)^ ^ {S^,S^)^ ^^^^ 
{h^hu) (S^Siy) - x'^ 

We now take the limit in the self-consistency condition {Tj). We keep Eg to derive a general formula, and we will 
specialize to various schemes later. We only use the fact that the diagonal part of Es is independent of K and equal to 
the cluster self-energy : Eg^^ = E^^. This is true for the schemes studied in this paragraph : PCDMFT, CDMFTa.nd 
DC A (For the pair scheme, see Appendix lD|) . Therefore, the dominant part of order U in Eg is diagonal and we have 

E5(X) = Ec;7 + 5E5(if) (32) 
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where ST,s is a matrix of order 0(1) in the l/U expansion. Denoting with a bar the normahzed sum of the reduced 
BriUouin zone, 



KeR.B.Z 



E 1 

KeR.B.Z 

we obtain [R is on the superlattice) : 

A.^x) = + I Jr(x) (hp) + O (^i^) (33a) 

i» = ti+JSi'-(5E^ (33b) 
Jr(-) ^ 7 {{ts + 6^s)Jts + SEs)^_^ {ts + 5J:s)Jts + SEs)^^) (33c) 



fl#0 

Note that t is purely off-diagonal. In the expression of J, (5Ss has to be expanded to order 0(1) only. Moreover, 
only the off-diagonal part (in site index) is important since the on site part is restricted to i? = 0. 5Y,s and A are 
determined by Eqs. (|29l I33|l and the relation between Eg and S*-'. 

At this stage, it is useful to distinguish two cases depending on the validity of the cancellation 

= E^ (34) 

In CDMFT and DCA , Es = E*^ is fsT-independent, therefore (PU holds and (5Es drops out of (|33c(l . i and J do not 
depend on x, although (/i^) does. In PCDMFT howeveT, (|34|) does not hold and we have to solve for (5Es and A (See 
below) to complete the computation of Gq. 

The second step of the computation is to take the large-?/ limit of the effective action H26|l using the value of Go 
()33|l. In that limit, we expect that the problem becomes classical, and more precisely for C7 —> oo: 

{cl^Ca„) iov V (35a) 
n-MT + '^Mi ^1 (35b) 

Indeed this can be shown explicitly using H27a|l and taking the J7 — > oo limit in the Fermi factors, after doing the 
summation over frequencies. However due to the frequency dependent nature of Go in all cases, we need to use the 
functional formalism and not the Hamiltonian formalism as in the derivation of the Ising limit on the lattice. This 
is presented in detail in Appendix IbI In the limit [/ — *■ oo with /3 — > oo, /3 = (3t/U fixed, the effective action reduces 
to the classical action for the cluster at temperature 1//3 and we obtain a classical cluster scheme for the Ising model 
with an interaction Jising inside the cluster and an additional mean field term Jb '■ 

Hoff = E Jisl,S,S, + E -^bS, {S.)h^^^ (36a) 

dx J/f'"(x) 
V(l-f-a;2)2 

where (/i, v) and (/i, v) correspond to a general and a nearest neighbors couple of sites respectively. One has to solve a 
classical Ising cluster, with self-consistent "boundary" condition represented by Jb, which generalize the usual Weiss 
field. Of course, if t does not depend on x, we find the same result as in the lattice Jj^h^g = J]f^y for nearest neighbors. 

Let us now specialize our computations to the three schemes presented in Section^ and compute the value of Jising 
and Jb for CDMFT, DCA and PCDMFT for square clusters of linear size L on a two-dimensional square lattice and 
for a general hopping ts where (5 is a lattice vector. 



J'b^I ^7T^ (36c) 



• CDMFT 
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In this case, T,s{K) = is K independent and ^5 = t, which leads to i{x) = t and : 

€-n, = (37a) 
= (-1)^ E J^-^+R (37b) 

The interaction inside the cluster is the same as in the lattice problem and the Jb term is of order 0(1) and is 
confined to the boundary of the cluster. The boundary term couples a spin to the average value of its "ghost" 
neighbor in the neighboring cells, this average value being computed in the cluster itself using the translation 
invariance on the superlattice. We added a — sign for odd cluster size, since in this case, the CDMFT has to 
be generalized like DMFT with two sublattices in order to capture antiferromagnetism : this is equivalent to 
reversing the sign of the "ghost" neighbor. In the large cluster limit L 00, the boundary terms play no role 
and we therefore recover the lattice Ising problem. Notice however that the one dimensional case is pathological 
since the two boundary terms communicate with each other resulting in a finite Tc in the limit of infinite size. 
This pathology disappears in higher dimensions and explains the results of reference (idj . 

DCA 

In (i = 2, we use a square cluster of linear size L — 2Lc, corresponding to points (Cf Appendix^. With 

the definition sinc(a;) = Y\i=i sin(xi)/xi : 



(38a) 

Ji'sL, (38b) 




7t{S - 6') 



5' =fi~i^[2Lc] 



Jb is the same for all links as required by translation invariance in the cluster. For the first neighbor hopping, 



denoting by J^^** the value of J^^** for i,j nearest neighbors, the formula reduce to : 



JiL. = ^J"-"'' = 1^2 - i| ) J^-" for = 1 (39a) 



>/Cg - (^Binc ^ j J"^" = ( 1 - ( sine — ) ) J"^" for > 1 (39b) 

Note that in the 2x2 cluster, the hopping is doubled since electrons can hop from one site to the neighbor either 
directly or using the cyclicity condition. In the large cluster limit Lc ^ 00, we recover the lattice problem: 
J^L^ ^ Jll'^'u and Jb = 0{l/Ll) mj. 

• PCDMFT 

This case is more complicated since Ss = S'^*' is now K dependent, hence ^ E^, except for diagonal 
elements. Therefore t and J now depend on x. In order to determine them, we have to solve Eqs. H9I29I33|I . 
This is done in Appendix |21 and we obtain : 



2 




{K) {K) 
' A{x) {h^,K) 



V(^)=Vir7ZT7prT (40a) 



t 



2 



Jr(^) = {-l)''A-\x) E (40b) 

P+sec 



where 5 is one of the basis vectors of the square lattice. 
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In d = 1 for first neighbor hopping, the last term reduces to 

The boundary terms are located as in CDMFT. In the large cluster limit L ^ oo, in Eq. H4U|) the boundary 
terms are subdominant in the denominator and the terms involving the averages of S cancel, restoring the 
correct Jismg on the lattice. The Jb term cancels since it is restricted to the boundary. 

• Pair scheme 

The computation of the classical limit of the pair scheme is a little more involved but similar. We show in 
Appendix ^ that if a magnetic solution exists, it satisfies the equation of the classical variation method : let's 
take two cluster problems, the first one with one site (denoted with an index (1)) in a field —zh, the second one 
with two sites (denoted with an index (2)), interacting antiferromagnetically with ji^^*' and in a field —{z — l)h, 
where z = 4 is the connectivity : 

= zhS'-^^ (41a) 

= _ _ 42)) ^ jLatt^(2)^(2) ^^^^^ 

where the field h is determined by the "nested consistency condition" 

< 5(1) >=< 5f ^ >= - < 5f ' > (42) 
Solving the classical equations, the critical temperature is given by : 

(z- I)tanh(/3JL^") = 1 (43) 

which gives f^/J^^^^ w 2.88. 



B. Discussion 



Let us recapitulate our results. We applied the cluster dynamical mean field theories to the Falikov Kimball model 
taking [/ — s- 00 with f3 — s- 00, P = (3t/U fixed. In this case the quantum model reduces to the Ising model and the 
cluster dynamical mean field approximation to classical cluster methods. CDMFT corresponds to a simple mean field 
theory generalized to clusters: each spin on the boundary is subjected to a mean field representing the interaction 
with the neighboring spins that do not belong to the cluster. The form of the mean field is the standard one: the 
antiferromagnetic coupling times the magnetization of the spins. Note that in principle the antiferromagnetic coupling 
is generated via quantum fluctuations and its value is approximation-dependent. In the case of CDMFT one gets the 
same coupling as the one obtained for the lattice. 

The classical limit of DC A is different mainly for two reasons: i) the mean field is not just on the boundary but 
it acts on all sites and is equal on sites belonging to the same sub-lattice. This is natural in DCA because it's an 
approximation that preserves the translation invariance (or the reduced translation invariance for the AF phase), 
in the propagator and in the Weiss field; ii) the form of the mean field is the standard one but the value of the 
antiferromagnetic coupling in the mean field and for spins inside the cluster are different from the lattice one (and 
between themselves). Of course in the limit of infinite cluster they reduce to the lattice value. 

PCDMFT is similar to CDMFT to the extent that the mean field acts only on the boundary. However, as for DCA 
, the value of the antiferromagnetic coupling in the mean field and for spins inside the cluster are different from the 
lattice one (and between themselves). Finally, the pair scheme reduces to the classical cluster variation method for 
two sites. 

The self-consistent equations corresponding to these classical cluster schemes can be solved analytically for small 
clusters and using a classical Monte Carlo method to solve the impurity problem for larger size. With the exception of 
the cluster variational method, the classical limits of these extensions of DMFT do not result in drastic improvements 
in the estimation of Tc for small cluster sizes. The value of Tc predicted by DMFT is the standard mean field one: 
4J (the connectivity of the lattice is four) which is far from the exact value that is 2.27J. The value obtained using 
the NCS scheme with two sites, i.e. the pair scheme, lead to a good improvement: Tc = 2.88 J. Instead the results for 
DCA and CDMFT even for larger sizes (4 by 4, 6 by 6 and 8 by 8 clusters) do not improve the estimate of Tc very 
much. In fact within error bars of 0.1 we get the Ti^^^'' = 3.2J, Tc^^^'^ = 3. J and yj*^^) = 2.85J (within the error 
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bars the estimate for Tc are the same for the two schemes). Note that the vahie obtained with a cluster of 16 sites 
with DC A and CI^MFT reaches the estimate obtained with a cluster of two sites using NCS! The results of PCDMFT 
are clearly worse than the ones of the two previous methods for the numerical value of Tc which is quite larger than 
what is obtained by CDMFT and DCA . This can be traced to the lack of cancellation of off-diagonal elements in 
(|34|l which implies that the cavity field A(icj„) in PCDMFT \s not proportional to the square of the hopping matrix 
element which is an undesirable feature. 



In this section, we present a general method to prove the causality of cluster approximations and apply it to various 
schemes, including PCDMFT. As mentioned in section there are two equivalent presentations for the equations 
of DMFT and its extensions : a first one using the Weiss function Go , and another one using the Luttinger-Ward 
functional Let's examine the causality question in both. 

• In the formulation with Gq, one has to prove that: i) If Gq is causal, Gc and S are causal : this is a automatic 
(as long as one uses a causal impurity solver), If E is causal. Go computed from the self consistency condition 
is causal. This is the difficult part. It was carried out explicitly for CDMFT in and for DCA in 0. 

• In the formulation with $, one has to prove that: i) Given a causal G, Yiskei{G) = 6^/6G is causal; ii) Given a 
causal E, the self consistency condition produces a causal G. Here ii) is obvious, because of the simple form of 
the bare propagator, i) is the difficult part and this section is devoted to a general method to prove it. 

To show that the self-energy is causal we have to prove that the retarded self-energy has a negative imaginary part 
for all w, i.e. ImT, ji {x, x' ,ui) = {Tir^x, x' ,uj) — T,*j^{x' , x,uj))/2i is negative. We use the Cutkovsky-t'Hooft-Veltman 
equation (also known in the literature as "cutting equations" ) which is extremely useful to discuss causality properties 
in terms of a diagrammatic expansion. It relates the imaginary part of the self-energy to a sum of cut diagrams. A 
standard derivation can be found in |2ll | and it is due to t'Hooft and Veltman. However, we present in Appendix IeI a 
simpler and self contained derivation, based on the Keldysh method. Contrary to the previous derivation it does not 
assume translation invariance, which is important to discuss cluster schemes that break translation invariance. Note 
however, that this method is limited to zero temperature formalism. 

This section is organized as follows : in paragraph II V Al we present the Cutkovsky-t'Hooft-Veltman equations and 
show how it can used to prove causality; in paragraph IIV Bl we apply it to various cluster schemes. 



The Cutkovsky-t'Hooft-Veltman equation gives the imaginary part of the retarded self-energy as a sum of cut 
diagrams. Let us consider the set of all perturbative connected diagrams {D} for the self-energy at zero temperature 
(for a given approximation) and all possible cuts of these diagrams into two connected parts L and R containing 
respectively the left and the right external point (if they are the same, we drop the diagram). There are n cut 
propagators going from left to right and n—1 from right to left, at frequency LOi and indices Xi. We denote by MR the 
mirror of R defined as the left diagram obtained from R by reversing all arrows. (Cf. Fig. Qfor an illustration). To 
each left diagram, we associate the value of corresponding T — diagram, using the T = Feynman propagator 
(G++ in Keldysh notation : see appendix IE)|. Note that for each given half left diagram one has to put in all 
the half left diagrams obtained permuting the cut lines in all the possible ways and then multiply by an appropriate 
symmetry factor discussed in the appendix (0. We then write DL{x,uj,{xi},{u!i}) = dLd{u!) + d*j^9{—Lo), where 9 
is the Heaviside function. We denote all cut propagators with an index i. To each we associate p{xi, x'^,uji)0{eiUJi) 
with ei — sgntj for a propagator going to the right, — sgnoj for a propagator going to the left. The summation over 
diagrams can be expressed as a summation over n and a sum over diagrams = {Ln, Rn) with a given number of 
cut propagators going to the right. The Cutkovsky-t'Hooft-Veltman equation is then: 



IV. CAUSALITY 



A. The Cutkovsky-t'Hooft-Veltman equation 




n>2 



D„ = (L„._R„) f=^i} 
{='} 



X 



X 



n 



Xi,x[,uJi)9{tiUJi) (44) 



i=l 
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FIG. 1: Definition of cut diagrams 



Remark: if the system is translation invariant then we can Fourier transform also with respect to space and get: 




(45) 

Note that there are no terms with n < 2, otherwise the original closed diagram would not be 2PI. 
Let us now see how to use (|44ll to prove causality of the self-energy. We want to show that: 

^.Y,wA^R{x,x',u)^^*j^ix',x,Lj))w:, <0 y{w,} (46) 
for all complex vector Wx- Using eq. (|44|l we find that we can write the previous quantity as a sum over cut diagrams: 

"R>0 {xj} 

Y[ p{xi,Xi,u!i)9{eiUJi) (47) 
1=1 

At this stage, we can use the causality of G, i.e. the positivity of p{x,x'). It implies the positivity of the tensor 
product nr=i^"" p{xi,x[,LUi) considered as a matrix of indices {xi, . . . ,a;„^+„^). Hence, positivity of — ImE would 
result if the summation over diagrams could be recast as a square. More precisely, a simple property that implies 
causality is that for any n : 

{(L„,i?„)} = {L„}x{i?„} (48) 

{L„} (resp. {Rn} ) is the set of left (resp. right) parts of diagrams with fixed n. Eq. H48|l says that all diagrams can 
be obtained once and only once using a left part and a right part, which is equivalent to two properties : 

• Closure property : With any left part and any right part (that have the same n), gluing them constructs a 
diagram present in the expansion. 

• Exact counting property : all diagram are obtained that way (by definition of left / right parts) but only once : 
there is no overcounting problem 

If the property H48|l holds, then using the notation 

Vn{{xi}\{uji},uj) = ^ w^i?L„(a;,Ci^,{xi},{a;i}) 
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we have : 

/riL+nB. 
X! yrii{^i}'d^'-}^^)yniWi}'A^i}^^) n P(2;<,a;-,wj)6'(e,w,) (49) 

Note that the sum over inside the integral can be interpreted as the product of the two vectors y and y* 

with a positive definite matrix. Thus the positivity of p imphes the negativity of ImE. 

The efficiency of this method is that one can prove causahty, just by examining the combinatoric structure of the 
diagrams present in the approximation, with no further computation. A few remarks are important at this stage : 
i) This method can prove the causahty of a scheme, but it can not prove that a scheme wih violate causality. In 
particular it proves the positivity of — ImS for all causal propagators G. In the following, we will refer to this as 
the strong causality property. It is a logical possibility that this property does not hold, while the actual solution of 
the scheme G is indeed causal for all parameters. However, we know no example of this, and the strong property is 
interesting since we use iterative algorithm to solve the self-consistent impurity problem, ii) Equation (|48|) is sufficient 
but not necessary : one could also have form squares in a more complicated manner. Hi) As a remark, we see that the 
full perturbation theory on the lattice satisfies the property (|48|l : the full diagrammatic expansion obviously leads 
to a causal self-energy. 

B. Applications 

1. CDMFT and DC A 

CDMFT can be obtained from the Baym-Kadanoff functional taking $ = ^{Ga.p:H..i{), where R is the cluster 
index and a is the internal cluster index. As a consequence Sq,/3;r.r = "^a^p = ^ ■ Since ^{Ga.p) is the sum of 

all the 2PI diagrams of a L'^ lattice, CDMFT does satisfy H48|l and therefore is a causal scheme. 

In DC A , one writes the self energy in terms of the complete series of diagrams and the delta function at the vertices 
are replaced by a coarse grained function. As discussed previously DCA is equivalent to replace the Kronecker delta 
Ski+k2.ka+ki with a Kronecker delta 5k^+k^,k^+k2 where fcf is the cluster momentum related to the hypercube containing 
ki. Formally kf — ki mod {t^/ . . . , tt/L). Hence, DCA satisfies (|48|) : the proof is the same that for the complete 
series of diagrams for the self energy. The fact that the k structure of the vertices is different does not change nothing 
at all. Hence, we see that DCA is automatically causal. 

2. PCDMFT 

The causality of PCDMFT is easily established using its real space formulation 0. The causality for CDMFT 
implies the causality of S*-^. We just have to show that S'''** is causal. For any complex vector Xa/j,, we have from ^ 

^ <^,x.vS!,^;V^,(i^) = J; E E 4:AS>.'Am''aZ', (50a) 
z^—^{5)^^^^e^^-V^\L (50b) 

The negativity of Yp implies that we have a sum of negative terms. This shows that causality of PCMDFT is 
guaranteed term by term in the S sum. 

Note that another way to reach the same conclusion is based on the fact that PCDMFT can be interpreted as 
the solution of a real lattice systems replacing the real U{ki, k2, k^, k^) with its PCDMFT counterpart. Thus, the 
statement of causality of PCDMFT is a particular case of the causality of the original problem with a general U. 

3. Non-Causality of the pair scheme 

The pair scheme, described in Section III CI has been originally introduced in 0, Q . It seems a very natural way 
to introduce systematically a k-dependence for the self-energy. Unfortunately, as we will show in the following, this 
scheme is not causal for two reasons that we shall discuss in detail. 
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First, let's apply the Cutkovsky rule. We see that the pair scheme does not satisfy the property (|48|l or more 
precisely the closure property. Indeed, when cutting a diagram for S, we only have two sites indices i and j. Let 
consider a cut diagram (LI, Rl) where we have i along the cut, and the similar diagram (L2, R2) where j is replaced 
by fc ^ A priori, we could glue the parts into {LI, i?2), but that diagram would involve i,j and k, and thus is not 
present in the diagrammatic expansion generated by the pair scheme. This is a first reason for lack of causality. We 
stress that this is not a proof, but a quick and strong indication that the scheme is not strongly causal. It's interesting 
to point out that this problem was already encountered in the case of disordered electron systems when people tried 
to generalize CPA to clusters. In that case Mills and collaborators 15] noticed that a generalization of CPA similar 
to the pair scheme has this type of non-closure problem. For this reason they introduced the traveling cluster method 
which is fully causal. We will comment afterward on the possibility to do the same for strongly correlated electrons. 

Let us now really show that the pair scheme is not causal. In fact the violation of the closure property is just 
an hint that there could be some problem with causality but it doesn't preclude the possibility that eventually the 
scheme is causal. In the following we will consider the solution of the pair scheme when J7 ^ for an arbitrary bare 
spectral density or hopping on the original lattice. We will show that a particular choice of the bare spectral density 
gives rise to a non causal self-energy thus proving in an explicit example the non causality of the scheme. 

To study the sign of the imaginary part of the self-energy in the limit C/ — > one can focus just on the second order 
term because the first order gives no contribution. Moreover one can replace in the diagram the full propagator with 
the bare or lattice one [Gbare = {i^n — t — the error leading to a term of an order larger than 0{U'^). If one 

applies the Cutkovsky equations to obtain the imaginary part of the self-energy for the second order diagrams of the 
pair scheme then one has to restrict the sum in the first line of Fig. |21to single sites {x,x) and links < x,x' >. The 
sum over the single sites {x, x) is just the same that one obtains for usual DMFT and is clearly positive. The link 
term <x,x' > with x,x' nearest neighbors equals (w is the (positive) self-energy frequency): 



The spectral density in the previous equations is just the bare spectral density obtained by tij. For the pair scheme 
one can find a spectral density matrix that, injected into the previous equations, does not lead to a causal self-energy. 
The causality of the spectral density matrix is equivalent to the positivity of each Fourier element {oj) > 0. However, 
even if pk{^) > 0, Px,x'{^) with x,x' nearest neighbors can be negative and as close as possible (in absolute value) 
to Px,x{^) for any oj. Consider for example the example in which pk{i^) is strongly peaked and has an important 
weight only near the points ±7r, . . . , ±7r. In this case px.xi^^) > 0, px.x'ii-^) with x, x' nearest neighbors is negative and 
Px,x'{^) = Jbz J^Pki^) exp{-ik{x - x')) ~ - /g^ j0)dPk{'^) = -px,x{^)- In t^is case the contribution of (EU, for 
example for = 1, is negative and approximatively two times larger in absolute value than the contribution from 
the x,x term. Since for each x,x term there are z {z is the connectivity) terms like H51|) the global sum is negative 
for z > 2, i.e. in any dimension large than one (if one takes a square or cubic or hypercubic lattice). Thus in this 
case, the self-energy is not causal. So we have found a second reason for non-causality: even if the closure property 
is verified and in the cutting-folding procedure no new diagrams are created (which is the case for the pair scheme at 
second order in U) this does not guaranteed that there is a way to write all the cut diagrams as a sum of squares. 
One can have an overcounting problem as we have just shown for the pair scheme. 

Let us now comment on the possibility to cure the pair scheme as Mills and collaborators did in the case of disordered 
electrons defining the traveling cluster scheme. What they did is to solve the violation of the closure property allowing 
an arbitrary self energy diagram to have only single site and nearest neighbor propagators. So in this case a self-energy 
diagram can connect two arbitrary points but just through a sequence of nearest neighbor propagators (that is why 
the scheme is called "traveling"). This clearly solve the first problem: the non-closure but in principle it doesn't 
guarantee that the second one also is cured. And indeed it is not for strongly correlated electrons. A simple way to 
see it is that even with a generalization of this type in the limit C/ — > the causality counter-example exhibited for 
the pair scheme is still valid. Instead the case of disordered electrons has a simpler diagrammatics (there are no loops 
in the self-energy diagrams due to the use of the replica method) and one can explicitly rewrite everything as a sum 
of squares. In particular the diagram discussed before is simply not there. 

A way to circumvent these two difficulties is write the self-energy as the sum of all diagram and then replace the 
original propagator with a simplified restricted version that is guaranteed to be causal. For example in order to have 
a causal generalization of the pair scheme, close in spirit to the traveling cluster, one can write: 




+ Wj:>wIp{x', X, LUi)p{x', X, UJ2) p{x' , X,UJ-UJ2- t^s)) 



(51) 
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FIG. 2: Example of cut diagram. MR means mirror image of R. 



where Go, Gi are fixed imposing the conditions G(0, . . . , 0, — G(0, . . . , 0, w) and G(7r, . . . , 7r,a;) = G(7r, . . . , tt, w), 
and G = (Gq ^ — S)~^ This guarantees that G is a causal propagator as long as G is causal. So when G is inserted in 
the complete series it will give rise to a causal self-energy. The problems with this scheme is that: (1) it seems that it is 
not <i>-derivable. The reason is that in all the diagrams there are only nearest-neighbor propagator. However since the 
diagrams connect in general two sites that are not nearest neighbor it is difficult to think how one can obtain it using 
a restricted set of diagrams for <I>. Note however that is also the case of the traveling cluster scheme for disordered 
electrons. (2) It is not clear what type of impurity solver, if any, can be used to sum the series of diagrams. In the 
case of the traveling cluster many simplifications arise that allow one to solve (2) but not (1). the smaller becomes 
the causality violation. Let us finally remark that the two mechanisms behind the non causality of the pair scheme 
apply also to the general NCS schemes even though we expect smaller violations of causality for larger clusters. 



V. CONCLUSION 



Cluster schemes are a promising method for studying strongly correlated electrons. Compared to the study of small 
finite size systems, they offer the advantage that the thermodynamical limit is taking from the outset, and hence one 
can hope for smaller finite size effects. As a prime example, we now know that even a one site cluster captures many 
non trivial features of the Mott transition. 

On the other hand, there is no unique generalization of the single site dynamical mean field approximation. Some 
generalizations view the DMFT and its resulting impurity (or multiple impurity models in the case of the NCS) as a 
trick for summming selected classes of diagrams. An alternative view is provided by the cavity construction and by 
functional methods, where an effective action is constructed to compute the correlation functions of selected coarsed 
grained or cluster degrees of freedom. Finally, one can view DMFT as a choice of boundary conditions for finite size 
studies of a small set of degrees of freedom, where the boundary conditions recognizes that those degrees of freedom 
are periodically repeated in a medium. While all those points of view, lead to the same set of DMFT equations in the 
single site case, they lead to different constructions for cluster extensions, and the strenghts and limitations of these 
extensions need to be explored. This paper contributed in this direction by clarifying two important aspects of these 
extensions, namely the conditions on the cluster scheme needed for the scheme to satisfy causality, and the reduction 
of the cluster schemes to spin cluster methods in the classical limit. 

First, we introduced PCDMFT which is a causal generalization of the scheme proposed by Lichtenstein and Kat- 
snelson |T^. Second, we provided a general way to prove causality of cluster scheme and showed that the pair scheme, 
which is the most natural extension in the sense that it is defined by a translational invariant restriction of the Baym 
Kadanoff functional, is not causal : the diagrammatic reason for the failure of this method was clarified. Third, we 
also showed formally, in the context of the Falikov Kimball model how the semiclassical limit of the different cluster 
methods reduce to classical spin cluster methods. Both DCA and CDMFT give comparable answers even though 
their Weiss fields have a very different form (in DCA the Weiss field is uniformly distributed inside the cluster while 
in CDMFT it is focused on the boundary) . On the other hand PCDMFT^ have an unphysical feature in the classical 
limit which results from a Weiss field which does not vanish as the square of the hopping matrix elements. The nested 
cluster schemes are clearly superior and provide a quantum generalization of the cluster variation method which gives 
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critical temperatures that converge rapidly with cluster size. We showed that the nested cluster schemes violates 
causality. Our proof suggests that the this violation is going to be more severe, when a substancial part of the lattice 
self energy has a range that exceeds the size of the maximal cluster used in the NCS. On the contrary when the range 
of the self energy is contained in the maximal cluster, the diagrams generated by the Cutkovsky procedure which are 
not contained in the NCS are small, which is consistent with the early numerical studies which indicated no violation 
of causality in the pair scheme at high temperature. This observation is in the same spirit as a recent study of a 
new (non-nested) scheme, the "fictive impurity model method" , which connected violation of causality with peaks in 
momentum space in the lattice self energy |22l |. 

An important lesson that can be learned from our results is that the choice of an optimal cluster scheme may 
depend of the property studied and the physical system at hand. Indeed our semiclassical analysis revealed PCDMFT 
is not accurate in insulating regimes, at least in the one that we focused on doing the semiclassical limit. Instead 
CDMFT performs rather well. We conjecture that CDMFT is accurate in insulat ing cases and PCDMFT in metallic 
cases. Indeed a recent numerical in one dimension work points in this direction (123). 

Another important result concerns the remarkable accuracy of the NCS scheme compared to the other ones. Even 
if its performance in a case with stronger quantum fluctuations remains an open question, we think that it would 
particulary interesting to try to cure the causality problem keeping the nesting idea inherent to the NCS scheme. 
Another route to follow is try to incorporate the cellular dynamical mean field ideas of defining impurity models in 
adaptive basis sets into a NCS scheme. In this way one could try to adapt the basis to the problem so that the 
resulting self energy is short range and the causality problem, even if not avoided, is sensibly reduced. 

Finally, while the paper focused on a few schemes the techniques developed here are quite general and may play an 
important role in selecting and optimizing cluster methods for specific applications. 

Acknowledgments 

We thank A. Lichtenstein, A. Georges, S. Florens, S. Biermann for useful discussions, and M. Jarrell for pointing 
to us Ref. This work was supported by NSF grant N° DMR -0096462. We acknowledge the warm hospitality of 
the KITP program on Reahstic Studies of Correlated Electrons, at UCSB (NSF grant PHY 99-07949) where part of 
this work was carried out, and the participants of this workshop for discussions. 

APPENDIX A: REAL SPACE FORMULATION OF DCA 

In this appendix, we present a real space formulation of DCA 8], which is useful for orders that break translation 
invariance. Let us consider first the translation invariant case. The hopping of the full lattice can be written using 
the reciprocal superlattice and the cluster indices (in real space) : 

V W = E e'^''+''^y^----"h{K + K,) (Al) 

Here (in d=:l), — < K < j-, Kc — n = 0, . . . , L — 1 is on the reciprocal lattice of the (finite) cluster. Using 
ts = i in 01, we get CDMFT. Using ts = with 

t^y^^K) ^j^H e^^^-^-^-'^MK + K,) (A2a) 

= e-'-^^^VWe'^^"'" (A2b) 

we get DCA . Indeed, t^y^ is cyclic in the cluster indices (by definition of AT^), and provided that we obtain a translation 
invariant (cyclic) solution we can diagonalize all matrices in the cluster to get Eq. 

Let us now consider the case where the solution of the real space formulation of DCA breaks the translation 
invariance in the cluster into a smaller invariance : the (big) cluster is divided into small clusters of linear size L, 
and we have translation invariance in the big cluster when the small clusters are thought as collapsed into one point. 
For example, for AF order in d = 1, L = 2 (See Figure^. We will then denote by 1 < a, /? < L the points in the small 
cluster, hy 1 < A, B < Lc the positions of the small clusters in the big clusters, and by A, B the positions of the big 
clusters on the full lattice expressed in units of the original lattice. Any point on the big cluster can be described with 
a couple {a, A) and any point on the lattice by a triplet (a, A, A). Similarly, we denote by — tt < /c < tt an element 
of the original lattice's Brillouin zone, by —ir/L < K < tt/L an element of the Brillouin Zone of the superlattice of 
the small clusters, and by —Tr/{LLc) < K < tt/{LLc) an element of the Brillouin Zone of the superlattice of the big 
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FIG. 3: Definition of big and small cluster in d = 1. L = 2, Lc — 3. 



clusters. Finally, the reciprocal lattice of the small cluster are Kc = 2|2-,n = 0, ...,L-l,the same for the big cluster 
considered as a cluster of Lc small clusters, in units of the original lattice, are Kc = XT^'P = 0? • ■ ■ i — 1- Hence 
the reciprocal lattice points of the big cluster considered as a cluster of LLc points are Kc + Kc- We solve the DC A 
in real space with the big cluster of size LLc ■ 

GaA,/3s(t^) = — ^ (A3) 

tilUu^.B^^^ = jzhr ^ ^^iKMi..+..-.,-.n)t^K + Kc + Kc) (A4) 



Li 



where we used Kc • xa = [27r] in the last equation. Because of the reduced invariance, we can diagonalise using Kc 
and obtain : 

G^p{Kc,u) = X: ^ _ ,c,e(^ ^ K^yKli„,) _ j:^^^^Kc,co) 
Using the unitary transformation : 

Gc.(3{Kc,uj) = e-''^-(^=-^^)G„0(K„c^) 

we obtain : 

Gc,a{Kc,iu)=y = ^ = (A6) 

Let us now concentrate on the AF order : L — 2 and for simplicity in dimension 1. In this case, Kc ~ 0, tt. We will 
now show that ljA6p is equivalent to the fc-space formulation presented in (denoted by a upperscript M). In the 
AF phase, correlations appears between k and k + tt and the self-consistency condition reads : 



K 



uj -t{K + Kc + Ti) 



-1 



where G and E are 2x2 matrices, which are non diagonal in the AF phase : 
Using the transformation 

CK + Ck+tt Ck — Ck+t: 

V2 '^"^^71~ 

Eq. I|A7|I is equivalent to (|A6I) . Therefore the two formulations of DCA are equivalent in the antiferromagnetic phase. 

The real space formulation of DCA share with CDMFT a property which is useful in practice : one doesn't have 
to anticipate the appearance of ordered phase, i.e. to adapt the cluster scheme for the order to be described. The 
order will show up automatically solving the real space DCA equations. Moreover, it can be more complex than an 
AF order, as long as the cluster is big enough to contain at least one unit cell. 
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APPENDIX B: LARGE-(7 LIMIT OF THE EFFECTIVE ACTION : PROOF OF (15^ 

Here, we present some details of the derivation of the large-JJ Umit of Eq. H26|l. in the hmit where U ^ oo with 
(3 ^ oo, P/U fixed. The effective action is 



ScS ^ - drdr' 

J Jo Jo 

rlJ 2f / 1 \ 

Go^uA'^n) = t—S,. - - jjJri^) {hp) + O [jj^j (B2) 

We are going to compute the expression of partial partition function Z[{?ip|}] (also denoted as Z[{S'p}]) as a function 

of the spins and and recover We first compute the expansion of lnZ[{np|}] at second order in t and first order 
in J. We obtain : 



s 

where denotes the sum over couples, ( )_g is an average of the Gaussian action at fixed 5. We see that 

Zo[{Sp}] = Zt=j^o[{Sp}] cx Jlcoshj'^Sp') (B4) 

p ^ ^ 

is independent of Sp (since Sp = ±1) and therefore it will be dropped in the following. Using the relations : 

(^Cp{iuj„)cl,{iujl^)'^ = ~Gfp,{iujn)S^„,u'^ (B5a) 



and the fact that in the AF phase i{—z) = t(z) since (Si) + {Sj) = 0, we can reduce the expression (|B3|) . extract the 
SiSj term and obtain finally : 

ln(Z[{5p}]) = -/3i7eff + 0(t', J') (B6) 

with i?eff given by Eqs. The additional term J^'' does not contribute to the Jising since it would give a 0(l/[/^) 

contribution. Moreover, since Jising ~ P/C/ we do not need to compute the XjU term of Yp . 

Finally, higher orders in the expansion do not contribute. Indeed, since .Z[{S'p}] is a gaussian action and since we 
compute In Z[{S'p}]/Zo, we only get connected diagrams, thus we have one sum over Matsubara frequencies. Each 
Matsubara sum gives a factor PU (in the large-t/ hmit J^uj replaced by ^ J d{xU)) and each G a factor . 
Therefore higher order terms are subdominant. 

APPENDIX C: CLASSICAL LIMIT OF PCDMFT 

In this Appendix, we present the solution of Eqs. H9I29I33|I for completeness. We restrict ourselves to nearest 
neighbor hopping. We start by solving at order 1 in the 1/U expansion. From (|29|l and (|33|l . we have ior ^ ^ v 
(A^p = 0): 

(KK) 



a:„ = <„-,. + *i:-- -AJ^^i^ (CD 

where a = 1, 2 is an index for the two (AF) solutions as in J^. Dividing the cluster into its two sublattices A and B 
with C = Ca^Cb and using 0, we have (for all <E C) : 

'^S'-p,p+, I E -^^''p+^ + l E ^^""L^s (C2) 

pGC^ pec^ 
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where /3 = a if ^ € Ca, (3 = a otherwise, since if p + 5 ^ C, the exponential depends on K and averages to zero. For 
fj,,^i + S£C, Eq. (|UT|l leads to : 

(^s> - " " i ('.?> £ (^?> 

for all /i, (5. The left hand side depends only on whether fi Cz Ca or not. Therefore, for all 6, we can define for fi G Ca 
and ^ e Cs : 



h?h'? 



(they are independent of and (). Eq. IjCSp then leads to : 

F^^ =ts + AsaF^^ + BscF^"" (C5a) 
F,"^ Aso^F^^ + BsaF^"" (C5b) 



which can be rewritten as 



(1 - ^^"-4 _ = i5 (C6a) 

(1 - Bs») ^ - Asa^F^^ = i5 (C6b) 

The determinant of the equations for F is given by J1q(1 ~ — Bsa), which does not vanish for a generic x. For 
\d\ > 1, = by assumption hence F'^^^^ — 0. In the classical limit, the two solutions are related by a spin-flip, so 
^ ^ does not depend on a for /i, v nearest neighbor, according to (|31|l . For \5\ — 1, the unique solution is obtained 



for equal Fs = f^"^/^ , which is given by : 



Fs = ^ (C7) 



1 



pfcC 



Using the square symmetry, one can pick up one 5 to compute the denominator. Hence, we obtain Eq. H40a(l and 
(|40c(l (the hopping is nearest neighbors). Note that the hopping is rcnormalised at order 1 (in 1/C/), but it is still 
restricted to nearest neighbors. In order to compute the J term, we now focus on t + JE^'^** which satisfies : 



It is a translation invariant quantity, which is also restricted to nearest neighbors. It has a formula analogous to (|A1|) 
for the hopping, and its value inside the cluster is given by its K average. Using (|33cf) to compute the mean field 
terms J//'^ , we see that the computation is similar than for CDMFT, with just a rcnormalisation of the hopping (note 
however that t + (JE'*^*' is not the effective hopping i defined as the order 1 term of A). In particular, the Jb terms 
have the same range as the original lattice hopping. Reporting l|C8p into H33c() leads to Eq. (|40a|l . which completes 
the derivation. 
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FIG. 4: Notations and reduced Brillouin zone for the pair scheme (shaded area) 

APPENDIX D: CLASSICAL LIMIT OF THE PAIR SCHEME 

The first task is to generalise the pair scheme in the presence of antiferromagnetic order. Denoting the two sublattices 
by A and B, we have the approximation for $ : 

%a^r ^ - ^) [Y. ^it^"] + ^^^'^^3^ + Y *2 [G„; , G,, , ] (Dl) 

\^eA jeB J <^> 

We denote by (1) and (2) the one site and the two site model respectively. We group the two one-site models into 
2x2 diagonal matrices, and denote with an index diag the diagonal part of a matrix. To write the self-consistency 
condition, we take as unit cell one site on A, one site on B. K is now a vector in the corresponding reduced Brillouin 
zone, presented on Figure^ The hopping and the lattice self-energy are given by : 

tAA{K)=tBB{K) = {) (D2a) 

tAB{K)=tC{K) (D2b) 

C{K) = 1 + e-*(^i-^"^) + e*(^i-^2) + e-*(^i+^2) (D2c) 

E:f/4(i^) = E'- = (1 - z)eW^ + zEi^i^ (D2d) 

T}X'^{K) = G(i^) (D2e) 
and the Green function by the usual formula (with 2x2 matrices). 



G={ioJn~t{K)-T}--^^{K,iuj^))-^ 

The scheme implies a consistency equations for the Green functions : 

G^^^ = G 
g(i) =D{G) 

where we denote by D the linear operator that restricts a matrix to its diagonal. 

First G04 is diagonal (with same proof as for the other schemes ) but here it is not trivial : 

{G^^}y\iuj^) = iuJn + ^['\iuJn) " Sf *(icj„) = ILU^ - Z (j^f^iuJ,,) ^ [l^n)) (D3) 

{Gf^y\iuj,,) = iuj,, + T.f(iuj^) - E'f *(ic^„) = zc^„ ~{z-\) (e^ (ic^„) - {iuj^S) (D4) 
Using previous notations, we have : 

Ai'^(x) 



Af)(x) 



(D5) 
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Let's turn now to the up electrons. Introducing the notations : 



A = D{A) 



and, using that, for d a diagonal matrix independant of if. Ad 
self-consistency conditions : 



— A d (for a — 1,2), we can expand the 



Ed) 



with the expansion (|33|l . to get 



2t 

c7 



1 

c72 



(a) ^ — (q)^ — (q) 



t 



tip 



pv 



i{K,x) =t{K) + T}''^\K,x) 
Clearly, p = i with the notation 2 = 1,1 = 2. Moreover, 



diag 
diag 



,(2) 
'diag 



.(2) 
'^diag 



and the J terms are given by 



J^''> ^\Y.~^pp{R)ipA~R) 



R^O 



which implies 



J(i) = z Jo 

/^) = (z-l)Jo 



Therefore the equations for A simplify into : 



A 



(1) 



.(2) 
'^diag 



- E 



(1) 



)+2{z- 1) Jo (4 



(2) 



Using that m*^^) = to^^^ we have (ft,- 



(2)\ _ UW 



at dominant order and therefore 



(a,z = 1,2) 



(D6) 
(D7) 

(D8a) 

(D8b) 
(D8c) 

(D9a) 
(D9b) 

(DlOa) 
(DlOb) 

(Dlla) 
(DUb) 

(D12a) 
(D12b) 

(D13) 



From HD5ID13I) . using the computation of the classical field from A presented in Appendix we obtain the classical 
variation method defined in H41I42|I . 

Strictly speaking, we only prove here that, in the large-C/ limit, if there is a magnetic solution, it obeys the CVM 
equations. We have not shown that such non zero solution exist. A priori, one could wish to push the semi-classical 
computation further and compute the classical field explicitely to check that it gives the values prescribed by the 
consistency of the classical equations. However, this is much harder to do than in previous schemes for the following 
reason : contrary to all other schemes studied in this paper, there is no cancellation between the diagonal part of 
j^iatt g^^^ y^c^ Therefore, we have to compute correction to (h) to order in order to compute A to order 1/U 

(the relation (|30a|l is a priori valid only at dominant order). Because f3 ^ U, we need A up to order to get the 

such correction (like we need A at order 1/U to get a classical field of order 1). All these difficulties are hidden in the 



.(2) 

''diag 



S(i) term in (IdT21 
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APPENDIX E: PROOF OF THE CUTKOVSKY-T'HOOFT-VELTMAN EQUATION 

In this appendix, we prove the Cutkovsky-t'Hooft-Veltman formula using the Keldysh method [2^. We first briefly 
present our conventions. 

1. Notations 

We denote with a + the upper contour (going from —00 to 00) and a — the lower contour. The definition of the 
Keldysh propagators (for fermions) is: 

-i <T<j>{x,t)(t>'i{x',t') > = G++{x,x',t-t') (EI) 

-i<f(j){x,t)(t)''{x',t')> = G—{x,x',t-t') (E2) 

i <(l)\x',t')(l){x,t) > = G+-{x,x',t-t') (E3) 

-i < 0(x,t)</.t(a;',t') > = G-+{x,x',t-t') (E4) 

After Fourier transformation (,f{uj) ~ J dte^'^* f (t)) , we have in particular the following useful relations : 

p{x, x' , uj) = ImGnix, x' , uj) (E5a) 

TT 

G+-{x,x',uj) = 2i7Tp{x,x',uj)e{~uj) (E5b) 

G-+{x,x',u;) = -2inp{x,x\u;)e{uj) (E5c) 

{G++ix,x',uj)y = -G—{x',x,uj) (E5d) 

To illustrate the diagrammatics, consider for example as an example a second order self-energy diagram in Fig. [S] It 





x.t 



x',t' 




+ 



FIG. 5: A self energy diagram and its Keldysh counterparts, with the Keldysh indices and the overall factor, 
corresponds to 

Gx,x' ,t-t'Gx,x' ,t-t' Gx',x,t'-t 

The Keldysh counterparts of the diagram are obtained by replacing each vertex by his corresponding Keldysh coun- 
terpart (See FigO. They corresponds to: 

z-^H — h — h /-^H — h — — — ^ f^^^ — ^ -X- — — — 

^x,x' ,t-t'^x,x' ,t-t'^x' ,x,t'-t ^x,x' ,t-t'^x,x' ,t-t'^x' ,x,t'-t ^x,x',t-'t'^x,x',t-t'^x',x,t'-t ' ^x,x' ,t-t'^x,x' ,t-t'^x' ,x,t'-t 



where G ' are the Keldysh propagators. 
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2. Proof 

The first point is to write the imaginary part of the zero-temperature retarded self energy in terms of the Keldysh 
components : 

Iml].(.) = (E6) 
It 

An important point is that this relation must hold for each diagram individually, in order to be compatible with any 
approximation considered as a diagram summation. Indeed at T = for w > 

2t 

and using l|E5blE5cp . it is sufhcient to prove that 

E++ + S"" + + 1!"+ = (E7) 

is verified for each diagram individually. A way to prove that, which is actually in Veltman's book|23|) consists in 
noticing that each diagram with a + at the largest time vertex cancels against the same diagram with a — at this 
largest time vertex. This happens because G++(t - t') = G^+(t - t') and G {t ~ t') = G^^{t - t') for t > t' so the 
only thing that changes between one diagram and the other is a minus sign (associated with the — vertex and not to 
the + vertex). This is called by Veltman the largest time equation. To unveil this relationship we use the following 
notation: circle a vertex if it corresponds to (f)~ and do not circle it if it corresponds to (j)'^ . Then the relationship 
with Veltman rules to obtain the largest time equation will become clear. 

To obtain the Keldysh expansion of ImS^, we can thus take the expansion of E (as defined with any perturbation 
theory) and put Keldysh indices on it (as explained in paragraph IE First, we will focus on w > 0. Then 

S_| (uj) = 0, we only need to compute E ^., i.e. the diagram with a + at the incoming vertex, and a — at the 

outgoing vertex. The crucial properties of zero-temperature diagrams is the following : due to the presence of 
functions in frequency for G"' and G ^ (Cf. paragraph IE 1 |l 1 . the frequency flows always from the + vertex to 
the — vertex. This implies that a diagram is zero if it contains a — vertex to which no external line is connected 
and that is surrounded by + vertices because of frequency conservation (the same is true if one replace + with — ). 
This cancellation generalizes to a connected set, or region, of — to which no external line is connected and that is 
surrounded by + vertices. Since we only have two external vertices, the non-vanishing diagrams are those with a line 
that cut propagators dividing them into a -I- region L on the left (connected to the incoming vertex) and a — region 
R on the right (connected to the outgoing vertex). This line is the cut of the diagram and the sum over Keldysh 
indices reduces to the sum over all possible cuts, (see Fig. P). In the L (resp. R) region, we have to use G"*"*" (resp. 
G— ). Since G—{x,x',u) = -(G++(a;',a;,a;))*, the R part gives D^/fl times (-l)number of internal lines^ ^^^^^ 
is the (left) diagram part obtained from R by inverting all the arrows. Finally, to each cut line going from right to 
left (left to right) is associated a G^ (G ^) propagator. Hence we obtain the third rule of the text : each cut lines 
going from left to right (right to left) is replaced by p{x,x' ,uj)9{uj) {p{x,x' ,uj)9{—ijj)). 

So there are two last things that remains to be clarifed in order to get the eq. H44|) : (1) that the symmetry factors 
are taking into account correctly and (2) that the phase between the half left and right diagram is the one leading 
to eq. 1)44(1 . Concerning the symmetry factors one has to be sure that all the counting is done correctly. In order to 
do this is useful to replace each half (or left) diagram with the sum of all half diagrams obtained from the original 
one permuting the cut lines in all the possible ways. The symmetry factor that one has to attach to them is exactly 
the one needed to get back the right symmetry factor for the original diagram by gluing together the half-right and 
the half-left diagrams. The "gluing" operations means attach to each i line (see fig. [2Jl of the half diagram the 
corresponding i line of the left diagrams. Indeed to get the symmetry factor of a given diagram one has to write all 
the topologically equivalent way to obtain the same diagram. This can be done starting from the left and from the 
right, i.e. writing all the different ways to get the same topologically equivalent left and right diagrams and then 
attach them in all the possible ways that give rise to the original diagram. Thus, the symmetry factor related to a 
half or left diagrams are just all the different ways (contractions) that can be used to create them. In this way the 
operation of "gluing" together a half and left diagram in which all the cut lines have been permuted produces n\ (n 
is the number of cut lines) times the original uncut diagrams. The 1/n! in eq. 1471 is there exactly to balance this 
redundant n\ term. 

Finally, let us focus on the phase between the right part of the diagram and -DJ/^ is 

a= (_l)'«(_i)''"«(-l)^=«(""+i)(-l)"'-p(_i)"(i)"-i (E8) 
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The first term comes from the sign in G~~ = — G"'"^*, with 1^ the number of internal hnes inside the right part 
R] the second comes from the — Keldysh vertices, with the number of vertices of type a in the right part; the 
third comes from the (— factor of each vertices at T = with Ua the number of outgoing Hnes in the vertices 
of type a : this factor change under conjugation; the fourth comes from the nioop broken loops in the cut diagrams 

(we restrict ourselves for fermions here); the fifth and sixth terms comes from the difference between G-\ and p 

(in the cut propagators going from left to right and right to left respectively). The number of cut loops is given by 
nioop — n ~ 1. Moreover, summing all lines ending at a vertex in the right part, we have : 

J2 ""-R"" = lR + n [2] (E9) 

a 

where [2] is the reduction modulo 2. Finally we get a — i which, combined with ljE6|) leads to the Eq. (|44|l . Similarly, 

for w < 0, we have ImE^ = /2i. The left part has — vertices, the right part + vertices. Using a similar analysis, 

we get a = —i, which leads to Eq. (|44|l . 
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